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ABSTRACT 

We present new photometric data of the transiting planet HAT-P-12b observed in 
2011. Our three transit curves are modelled using the JKTEBOP code and adopting 
the quadratic limb-darkening law. Including our measurements, 18 transit times span- 
ning about 4.2 yr were used to determine the improved ephemeris with a transit epoch 
of 2,454,187.85560±0.00011 BJD and an orbital period of 3.21305961±0.00000035 d. 
The physical properties of the star-planet system are computed using empirical calibra- 
tions from eclipsing binary stars and stellar evolutionary models, combined with both 
our transit parameters and previously-known spectroscopic results. We found that the 
absolute dimensions of the host star are Ma=0.73±0.02 M©, i?A=0.70±0.01 R©, log 
#A=4.61±0.02, p A =2.10±0.09 p Q , and L A =0.21±0.01 L . The planetary companion 
has M b =0.21±0.01 M Jup , i? b =0.94±0.01 R Jup , log 5b =2.77±0.02, /> b =0.24±0.01 /> Jup , 
and T oq =960±14 K. Our results agree well with standard models of irradiated gas giants 
with a core mass of 11.3 M©. 

Subject headings: planetary systems — stars: individual (HAT-P-12) — techniques: 
photometric 



1. INTRODUCTION 

Transiting exoplanets are interesting objects to study because the coupling of radial velocity 
and photometric measurements allows a determination of stellar and planetary parameters and 
thus give us an important constraint on fundamental models of planet formation and evolution. 
For that purpose, we need to precisely measure the physical properties of both the stars and 
planets in these systems from a detailed analysis of transit curves obtained with high-precision 
photometry. Basically it is possible to regard the star-planet systems as detached eclipsing binaries 
with extremely low mass ratios (< 0.001). Therefore, the methods and tools that are used to model 
the binary stars could be applied to analyze the light curves of transiting exoplanets and then to 
determine their physical parameters (Southworth 2008, 2009; Lee et al. 2011). 

The transiting nature of HAT-P-12b was detected with the HAT-5 telescope of Hungarian- 
made Automated Telescope Network (HATNet; Bakos et al. 2004) in 2006. Hartman et al. (2009) 
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conducted follow-up photometry of four transits together with spectroscopic observations and re- 
ported that the transit features come from a Saturn-mass planet with a mass of 0.21±0.01 Mj up 
and a radius of 0.961q'q2 Rjup in a 3.2-d circular orbit. Its host star is a K4 dwarf GSC 3033- 
706 (2MASS J13573347+4329367; V=+12M) with M A =0.73±0.02 M , i? A =0.70±°;° 2 R , and 
[Fe/H]=— 0.29±0.05. Since HAT-P-12b is one of the lowest-density planets orbiting metal-poor 
host stars, the physical properties of the system are important for irradiation models. In this 
work, we report and analyze three new high-precision transits of HAT-P-12 and refine the physical 
parameters of the transiting planetary system. 

2. NEW OBSERVATIONS 

We obtained new photometric data from observing three transits of HAT-P-12b. Observations 
were carried out between 2011 March and May, using an ARC 4K CCD camera and a Cousins R 
filter attached to the 1.0-m reflector at the Mt. Lemmon Optical Astronomy Observatory (LOAO) 
in Arizona, USA. The telescope was significantly defocused, because it is expected to minimise 
random and flat fielding errors (see, e.g., Southworth et al. 2009). The e2v CCD chip has 4096x4096 
pixels and a pixel size of 15 //m. The image field-of-view is ~ 28x28 arcmin 2 at the f/7.5 Cassegrain 
focus of the telescope. A summary of the observations is given in Table 1, where we present 
observing interval, filter, binning mode, exposure time, numbers of observed points, and weather 
condition. During the first run, a gap before the transit ingress was caused by technical problems. 
With the customary IRAF package, we processed the CCD frames to correct for bias level, dark 
noise, and pixel-to-pixel inhomogeneities of quantum efficiency (flat field correction). We applied 
simple aperture photometry to obtain instrumental magnitudes. 

For each transit event we constructed an artificial comparison star by monitoring field stars 
imaged on the chip. Following the method described in Lee et al. (2011), we selected and combined 
by a weighted average about 10 candidate non-variable comparison stars. Then, the differential 
magnitudes from the artificial reference were normalized by fitting a linear function to the out-of- 
transit data to remove time- varying atmospheric effects (detrending). Resultant transit curves are 
plotted in Figure 1 and listed in Table 2, where times are Barycentric Julian Dates (BJD) in the 
Barycentric Dynamical Time system (Eastman et al. 2010). 

3. LIGHT-CURVE ANALYSIS AND TRANSIT TIMES 

To determine the planetary and orbital parameters, three LOAO transits of HAT-P-12 were 
analyzed simultaneously in a manner almost identical to that for the transiting planetary system 
TrES-3 (Lee et al. 2011) using the JKTEBOP code (Southworth et al. 2004a,b), which is a 
code for modelling the light curves of detached eclipsing binary stars using biaxial spheroids. The 
main parameters of the model are the orbital period (P), the ephemeris epoch (To), the fractional 
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radii of star (r\=R\/a, where a is the orbital semi-major axis) and planet (r\ ) =R\ ) /a), the orbital 
inclination (i), and the limb-darkening coefficients (LDCs). Actually, ta and are incorporated 
as their sum (r^+r^) and ratio (k=r\ 3 /rA)- Throughout this paper, we refer to star and planet as 
the subscripts 'A' and 'b', respectively. 

In the transit analysis, we used rA+rh and k as the fitting parameters, because these parameters 
are more weakly correlated than between r\ and (Southworth 2008). Initial quadratic LDCs 
were taken from the tables of Claret (2000), using the atmospheric parameters of Ta=4650±60 K, 
log # A =4.61±0.01, and [Fe/H]=-0.29±0.05 (Hartman et al. 2009). Final results are obtained from 
fitting the linear LDC (ua) but fixing the non-linear LDC (va), because the two LDCs suffer from 
strong correlations between them (Southworth 2008; Johnson et al. 2008). These are summarized 
in Table 3, together with the stellar density (pa) and the planetary surface gravity (g\>) and zero- 
albedo equilibrium temperature defined as T eq =TA\/ ta/2 (Southworth 2010). The lower three 
values were directly calculated from those transit parameters and the stellar velocity amplitude 
(Xa=35.8±1.9 m s" 1 ) of Hartman et al. (2009). In order to assess the uncertainties of the fitted 
parameters, we ran 10000 Monte Carlo simulations and a residual permutation algorithm (Jenkins 
et al. 2002) implemented into JKTEBOP, respectively. On these occasions, the non-linear LDC 
was perturbed by ±0.1 around the fixed value. The error estimates presented in Table 3 are the 
lcr-values adopted from the larger of the two results. Figure 2 displays the light curves with our 
best-fitting model and residuals. 

For ephemeris computations, the minimum times and their uncertainties for our three transits 
were determined with the JKTEBOP code and the model parameters of Table 3. In addition 
to these, four transit times were newly determined by us from the individual measurements of 
Hartman et al. (2009). The times reported in the discovery paper are not actual events, because 
they were calculated by assuming strictly periodic transits and by adjusting both centers of the 
first transit and the last transit covering all observations in transit modeling. The results are given 
in Table 4, together with 11 epochs taken from Exoplanet Transit Databasqj (ETD; Poddany et al. 
2010). The ETD timings were transformed from HJD based on UTC into TDB-based BJD using 
the online applets^ developed by Eastman et al. (2010). By including these literature-based transit 
timings into a linear least-squares fit, we computed the orbital ephemeris of HAT-P-12 to be 

C tr = BJD 2, 454, 187.85560(11) + 3.21305961(35)^, (1) 

where E is the number of orbital cycles elapsed from the reference epoch and the parenthesized 
numbers are the lcr-error values for the last digit of each term of the ephemeris. Individual mid- 
transit timing errors were used to obtain the corresponding weight factors. The period is somewhat 
shorter than those of Hartman et al. (P=3.2130598±0. 0000021 d) and from our transit analysis in 
Table 3. 



1 http://var2. astro.cz/ETD/ 

2 http: / /astroutils. astronomy.ohio-state.edu/time/ 
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The observed (O) minus calculated (C tr ) transit times from equation (1) are given in the 
fourth column of Table 4 and plotted in Figure 3. The timing residuals seems to indicate a short- 
term oscillation as is shown in the figure. This might be giving a hint about a period change, for 
example, due to the presence of a second planet in the HAT-P-12 system. We examined whether the 
residuals represent periodic variability, but no detectable periodicity was found. As in the case of 
close binaries with deep convective envelopes (Lee et al. 2009), the apparent transit timing variation 
in Figure 3 may be produced by the sporadic asymmetries of transit curves due to starspots on 
stellar surface caused by magnetic activity. In reality, previous measurements (especially from 
ETD) might have red noise contamination and their errors would then be underestimated (Carter 
& Winn 2009). 



4. RESULTS AND DISCUSSION 

We can compute the absolute demensions of transiting planetary systems by combining pho- 
tometric and spectroscopic results. Currently, there are two main methods to obtain stellar and 
planetary parameters from the observed quantities. The first method is to use the empirical rela- 
tions from eclipsing binary stars and the second to apply the isochrones from stellar evolutionary 
models. For this procedure, we used the planet velocity amplitude as a key parameter governing 
the solution process to find the best match between the observations and predictions, which is the 
same approach as in Southworth (2009) and Lee et al. (2011). 

First of all, we calculated the physical properties of the HAT-P-12 system using the new 
calibrations of stellar masses and radii expressed as Ta, log pa, and [Fe/H] (Enoch et al. 2010, 
hereafter ECPH), which are originally defined by Torres et al. (2010) from 95 well-studied eclipsing 
binaries and replaced log qa with log pa- The process is to look for K^, satisfying simultaneously 
the stellar mass and radius from the two relations. The calibrated mass and radius of Ma,ecph= 
0.765±0.012 M and -Ra,ecph= 0.686±0.012 R© were used to minimize the \ 2 expressed as: 



x 2 



-^A,ECPH - MA jP red 



0"M 



A, ECPH 



2 



+ 



-Ra.ECPH - -RA.pred 



A, ECPH 



(2) 



where MA iPre d and i?A,pred are the relation-predicted mass and radius from the observations and 
the values. The results are given in the second column of Table 5. The quantity denotes the 
Safronov (1972) number. 

By considering both the metallicities allowed by the observational errors in [Fe/H] and the 
ages for each metallicity, we use three different sets of the stellar evolutionary models: Yonsei-Yale 
(Demarque et al. 2004), Padova (Girardi et al. 2000), and Baraffe et al. (1998, hereafter BCAH). 
This method also aims to find the velocity amplitude for which the measured values [Fe/H], 
ta, and Ta are best fitted to the radius i?A,pred and temperature TA iPre d predicted from the model 
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isochrones for each metalicity. This consists of calculating the \ 2 fitting statistic, 



x 2 



[Fe/H] - [Fe/H] model 



a 



[Fe/H] 



2 



+ 



r A - (-RA.predA 



>r A 



2 



+ 



T A -T A 

,pred 

<JT A 



2 



(3) 



where <7[F e /H]> a r A , and ctt a are the uncertainties corresponding to the measurements. Our process 
obtained a best-fit model when the BCAH isochrones with [Fe/H]=— 0.25 were used; both the 
Yonsei-Yale and Padova isochrones cannot constrain the age of the planetary system and favor the 
metal-rich stellar models. 

The isochrones from the BCAH models are plotted in Figure 4 along with the position of 
HAT-P-12A and the results are listed in the third column of Table 5. The mass and radius from 
BCAH are somewhat smaller than those derived from the empirical calibration of eclipsing binary 
stars, although the physical parameters between the two methods agree with each other within their 
uncertainties. Similar situations are found for WASP-21 (Bouchy et a. 2010; Barros et al. 2011), 
WASP-37 (Simpson et a. 2011), and WASP-39 (Faedi et al. 2011). As is the case for HAT-P-12, 
the three planet host stars are metal-poor with metalicities [Fe/H] of — 0.46±0.11, — 0.40±0.12, and 
— 0.12±0.10, respectively. On the other hand, 11 eclipsing binaries sampled by Torres et al. (2010) 
are low-metalicity systems ([Fe/H]<0.0). Among these only one component star (V636 CenB) is 
smaller than 1 M Q . From this it follows that the physical properties from the BCAH model for 
metal-poor low-mass stars seem to be more reliable than those from the empirical calibrations. We 
chose the BCAH solutions as our final results of HAT-P-12. 

The location of HAT-P-12b in the mass-radius diagram is shown in Figure 5, together with 
ten known Saturn-mass transiting exoplanets with masses in the range 0.15 Mj up < M < 0.4 Mj up , 
wherein Kepler-16b is a circumbinary transiting planet on a nearly circular 229-d orbit around its 
two parent stars (Doyle et al. 2011). In the same figure, we show constant density contours for 0.1, 
0.25, 0.5 and 1.0 pj up - The physical properties obtained in this study indicate that HAT-P-12b is 
a low-density sub-Saturn mass planet with a mass of 0.21 Mj up , a radius of 0.94 Rj up , and a mean 
density of /?b=0.24 pj up . The results are most similar to those (Mb=0.20 Mj up , R^=1.00 Rj up , 
Pb=0.19 Ajup) of the transiting exoplanet HAT-P-18b (Hartman et al. 2011). However, HAT- 
P-12 is younger and more metal-poor than HAT-P-18 with an age of 12.4 Gyr and a metallicity 
of [Fe/H]=+0.10. Furthermore, a classification of transiting close-in planets has recently been 
suggested by Hansen & Barman (2007). With the derived equilibrium temperature of T eq =960 K 
and Safronov number of 0.023 would classify HAT-P-12b to be of class II. However, it was recently 
shown by Southworth (2010) that a correlation between T eq and O appears to have little statistical 
significance. Therefore this classification seems to be of little importance. In addition, future 
observations of secondary transits (or occultations) of HAT-P-12b might prove to be difficult due 
to the planet's low equilibrium temperature. 

We compared our mass and radius to the predicted values from theoretical models of Fortney 
et al. (2007) with various core masses. For these, their models were interpolated to the age of 
3.2 Gyr and the solar equivalent semi-major axis of 0.0843±0.0029 AU calculated from ay^A and 
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were plotted in Figure 5 as solid curves. We conclude that HAT-P-12b is a H/He-dominated gas 
giant planet with a core mass of 11.3^2i and is moderately irradiated by its low-metallicity 
host star. Of the 10 circumstellar transiting planets, eight exoplanets including HAT-P-12b follow a 
suggestive correlation between the inferred core mass and host star's metallicity (Guillot et al. 2006; 
Burrows et al. 2007), while the recently discovered planets HAT-P-18b and HAT-P-19b (Hartman 
et al. 2011) with negligible core masses but super-solar metallicities disagree with the prediction. 
Further discoveries will help to identify and understand the possible correlation between planetary 
parameters such as metallicity, core mass, radius, and equilibrium temperature. 

The authors wish to thank the staffs of LOAO for assistance with our observations. We 
appreciate the careful reading and valuable comments of the anonymous referee. We have used the 
Simbad Database maintained at CDS and the Exoplanet Transit Database in this research. This 
work was supported by the KASI (Korea Astronomy and Space Science Institute) grant 2012-1- 
410-02. 
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Fig. 1. — Transit light curves of HAT-P-12 observed between 2011 March and May. 
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Fig. 2. — Phased light curves of HAT-P-12. The continuous curves represent the solutions obtained 
with the best-fit parameters listed in Table 3. The residuals from the fit are offset from zero and 
plotted at the bottom in the same order as the transit curves. 
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Fig. 3. — (O — Ct r ) residuals for HAT-P-12b constructed with the transit ephemeris (1) newly 
derived in this work. The minimum times calculated from our transits are plotted using the filled 
squares. 
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Fig. 4. — Isochrones from the BCAH models for log ages = 8.0 (bottom), 8.3, 8.7, 9.0, 9.3, 9.5, 9.6, 
9.7, 9.8, and 9.9 yr (top). The observed values of Ta and rx{=R^/a) for HAT-P-12A are shown 
together with their error bars. 
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0.1 0.2 0.3 0.4 



Planet mass (M Jup ) 

Fig. 5. — Mass-radius diagram of Saturn-mass transiting planets. W=WASP, H=HAT, K=Kepler, 
C=CoRoT, and symbols S, N, and U mark the location of Saturn, Nepturn, and Uranus. Dotted 
lines are density contours for 0.1, 0.25, 0.5, and 1.0 /9j up . Overlaid are Fortney et al. (2007) models 
interpolated to the solar equivalent semi-major axis of 0.0843 AU and the age of 3.2 Gyr for HAT- 
P-12b with core masses of 0, 10, and 25 M®. Data are taken from |http:/ /exoplanet.eu/] except for 
our result for HAT-P-12b. 
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Table 1. Observing log of HAT-P-12. 



Transit 


UT Date 


Observing Interval 


Filter 


Binning 


Exposure 


N ohs 


Weather 






(BJD+2,455,000) 




Mode 


Time (s) 




Condition 


1 


2011 03 29 


649.67-649.95 


Rc 


2x2 


50-70 


224 


clear 


2 


2011 04 14 


665.74-665.97 


Rc 


2x2 


50-70 


273 


partly cloudy 


3 


2011 05 13 


694.65-694.93 


Rc 


2x2 


40-60 


314 


clear 



Table 2. i? c -band photometry of HAT-P-12 



BJD 


Diff. 


Mag 


CMag 


Relative Flux 


2455649.67120 


+0. 


.0005 


0.0021 


0.9995 


2455649.67207 


-0. 


.0043 


0.0021 


1.0039 


2455649.67293 


-0. 


.0004 


0.0021 


1.0003 


2455649.67380 


+0. 


.0009 


0.0021 


0.9991 


2455649.67467 


+0. 


.0019 


0.0021 


0.9982 


2455649.67554 


+0. 


.0011 


0.0021 


0.9990 


2455649.67641 


+0. 


.0011 


0.0021 


0.9990 


2455649.67729 


-0. 


.0038 


0.0021 


1.0035 


2455649.67815 


-0 


.0001 


0.0021 


1.0001 


2455649.67902 


+0. 


.0001 


0.0021 


0.9999 



Note. - This table is available in its entirety in a 
machine-readable form in the online journal. A portion 
is shown here for guidance regarding its form and con- 
tent. 
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Table 3. Transit parameters of HAT-P-12. 



Parameter 


Value 




0.0969±0.0012 


k (=r b /r A ) 


0.1370±0.0019 


i (deg) 


89.915±0.098 


U A 


0.739±0.069 


VA 


0.177 perturbed 


T (BJD) a 


649.79751±0.00036 


P(d) 


3.213089±0.000037 


TA 


0.0852±0.0012 




0.01168±0.00005 


a (mmag) 6 


2.1269 


^rcd 


0.9920 


PA (P ) 


2.100±0.089 




0.257±0.014 


T cq (K) 


960±14 



BJD 2,455,000 is suppressed, 
'rms scatter of residuals. 
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Table 4. Observed transit times of HAT-P-12. 



BJD 

(2,450,000+) 


Uncertainty 


E 





-Ctr 


References* 


4,187.85655 b 


±0.00020 





0. 


.00095 


Hartman et al. (2009) 


4,216.77265 b 


±0.00014 


9 


-0. 


.00049 


Hartman et al. (2009) 


4,869.02397 b 


±0.00017 


212 


-0. 


.00027 


Hartman et al. (2009) 


A RQ7 Q/IOO^b 
4,83 1 .y^ZZiO 


_i_n nnn9/i 

ztU.UUUZ^i 


991 


u. 


nnn/is 

UUU4o 


jrLaxTinaii ex ell. t^zuuyj 


4,952.56398 


±0.00080 


238 


0. 


00019 


Gregorio (AXA) 


4,965.41639 


±0.00046 


242 


0. 


.00036 


Kucakova (TRESCA) 


4,965.41748 


±0.00090 


242 


0. 


.00145 


Ayoimanitis (AXA) 


4,984.69368 


±0.00060 


248 


-0 


.00070 


Gary (AXA) 


5,312.42673 


±0.00032 


350 





.00027 


Vilagi & Gajdos (TRESCA) 


5,630.51896 


±0.00049 


449 


-0 


.00040 


Ivanov & Sokov (TRESCA) 


5,646.58477 


±0.00059 


454 


0. 


.00011 


Nicolas (TRESCA) 


5,646.58486 


±0.00040 


454 


0. 


.00020 


Gajdos & Vilagi (TRESCA) 


5,649.79769 


±0.00020 


455 


-0 


.00003 


This paper (LOAO) 


5,659.43563 


±0.00038 


458 


-0. 


.00127 


Ruiz (TRESCA) 


5,665.86234 


±0.00031 


460 


-0 


.00068 


This paper (LOAO) 


5,675.49947 


±0.00064 


463 


-0. 


.00273 


Zinoviev & Sokov (TRESCA) 


5,694.78089 


±0.00024 


469 


0. 


.00033 


This paper (LOAO) 


5,704.42185 


±0.00038 


472 


0. 


00211 


Salas (TRESCA) 



a AXA (Amateur eXoplanet Archieve), TRESCA (TRansiting ExoplanetS and 
CAndidates). 

b Newly determined by us from the individual measurements. 
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Table 5. Physical properties of the HAT-P-12 system. 



Parameter 


ECPH Model 


BCAH Model 


K h (km s- 1 ) 


131.5±3.0 


129.7±1.5 


M A (M ) 


0.757±0.038 


0.727±0.019 


Ra (Re) 


0.711±0.019 


0.702±0.013 


log g A (cgs) 


4.613±0.032 


4.607±0.020 


La (L ) 


0.212±0.016 


0.206±0.013 


M boliA (mag) 


6.434±0.080 


6.463±0.069 


M b (M Jup ) 


0.216±0.015 


0.210±0.012 


Rb (Rjup) 


0.949±0.022 


0.936±0.012 


Pb (PJup) 


0.236±0.016 


0.240±0.012 


e 


0.0233±0.0020 


0.0236±0.0015 


a (AU) 


0.03887±0.00088 


0.03829±0.00046 


Age (Gyr) 




3.2±3.8 



